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ics, we derive the complete result for the two-loop correction to order {a/n)^ (Z a)"^ for the g factor of an 
electron bound in an nS state of a hydrogenlike ion. The results obtained significantly improve the accuracy of 
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I. INTRODUCTION 

The g factor of a bound electron is the coupling constant 
of the spin to an external, homogeneous magnetic field. In 
natural units /i = c = eo = 1, it is defined by the relation 



2 m \ / 2 
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where 5E is the energy shift of the electron due to the interac- 
tion with the magnetic field B, m is the mass of the electron, 
and e is the physical electron charge (e < 0). The matrix a ■ B 
contains the Pauli spin matrices a and has eigenvalues ±\B\. 

Studies of the free-electron g factor play an important role 
in modern physics. Together with the discovery of the Lamb 
shift in hydrogen, the observation of the electron magnetic 
moment anomaly led to the development of quantum electro- 
dynamics (QED). After decades of intensive theoretical and 
experimental studies, the free-electron g factor provides one 
of the most accurate and stringent tests of QED 1 1 ]. With the 
increased experimental and theoretical precision, it presently 
yields the most accurate determination of the fine-structure 
constant a 101 • 

It has not been until recently that investigations of the 
^jownii-electron g factor came into prominence. As was 
demonstrated in Ref. JsJ], the theoretical value of the bound- 
electron g factor can be used for the determination of the mass 
of the electron when combined with an experimental value for 
the ratio of the electronic Larmor precession frequency t^L and 
the cyclotron frequency of the ion in the trap ujc. 
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where q is the charge of the ion and mjon is its mass. The ac- 
curacy of the best experimental results for light hydrogenlike 
ions 01 is already below the 1 part per billion level and is 
likely to be improved in the future. According to the recent 
adjustment of fundamental constants these measurements 
provide the most accurate method for the determination of the 
electron mass. 



In order to match the experimental precision achieved, vari- 
ous binding and QED corrections to the bound-electron g fac- 
tor have to be calculated. It has been found long ago |6| that 
in a relativistic (Dirac) theory, the g factor of a bound elec- 
tron differs from the value g = 2 due to the so-called binding 
corrections. For an nS state, they are given by 
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where E is the Dirac energy. Other corrections to the bound- 
electron g factor arise from the QED theory. They were the 
subject of extensive theoretical investigations during the last 
decade. Accurate calculations of the one-loop self-energy 
iItI 11 I^JT3Jll|l^ vacuum-polarization |8, 9,"l3Cl, nuclear- 
recoil lll4lll5lfl6 1. and nuclear-polarizability flTl] corrections 
have been carried out. Detailed g factor investigations have 
been performed also for other systems that could be of experi- 
mental interest in the near future, in particular, for Li-hke ions 
and hydrogenlike ions with a non-zero nuclear spin fl^ . 

The subject of this work is the two-loop QED correction, 
which is presently the main source of the uncertainty of the- 
oretical predictions for the g factor of hydrogenlike ions. We 
present a complete calculation of this correction up to the or- 
der of (oLjiif' {Z otf. This two-loop correction has already 
been addressed to in our former work 1I20I1 . where an incom- 
plete calculation using a photon-mass regularization was pre- 
sented and an estimate for the total contribution up to the or- 
der (a/7r)^ (■2' ck)"* was obtained. The present computational 
method is based on the dimensionally-regularized nonrela- 
tivistic quantum electrodynamics (NRQED), which is a rel- 
atively new and very powerful approach for the calculation of 
higher-order relativistic and radiative effects. It has already 
been successfully applied to several challenging problems, 
e.g., to the calculation of the positronium hyperfine split- 
ting ll2lll and the ground-state energy of the helium atom I22ll . 
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II. DIMENSIONALLY REGULARIZED NRQED 

As is customary in dimensionally regularized QED, we here 
assume that the dimension of the space-time is £> = 4 — 2 e, 
and that of the space d = 3 — 2 e. The parameter e is consid- 
ered as small, but only on the level of matrix elements, where 
an analytic continuation to a noninteger spatial dimension is 
allowed. 

Let us briefly discuss the extension of the basic formulas 
of NRQED to the case of an arbitrary number of dimensions. 
The momentum-space representation of the photon propaga- 
tor preserves its form, namely g^ij/k"^. The Coulomb interac- 
tion is lEll 
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where the latter representation provides an implicit definition 
of Zg, and we have used the formula for the surface area of a 
d-dimensional unit sphere 
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The nonrelativistic Hamiltonian of the hydrogenic system is 



(5) 



Hn = 



p 
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The operator p ^ is well defined in any integer dimension. 
If we restrict our consideration to the spherically symmetric 
states, can be continued to an arbitrary real dimension by 
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In the following, we will not need the explicit (unknown) form 
of the solution of the Schrodinger equation in d dimensions. 
It will be sufficient to use instead its scaling properties, which 
we obtain by introducing the dimensionless radial variable p 
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In atomic units, i.e. expressed as a function of the dimension- 
less p, the Schrodinger Hamiltonian takes the form 
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We now turn to relativistic corrections to the Schrodinger 
Hamiltonian in an arbitrary number of dimensions. These cor- 
rections can be obtained from the Dirac Hamiltonian by the 
Foldy-Wouthuysen transformation. In order to incorporate a 
part of radiative effects right from the beginning, we use an 
effective Dirac Hamiltonian modified by the electromagnetic 



form factors Fi and F2 (see, e.g., Ch. 7 of ||23l]) 
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where 
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We use three-dimensional notations here, namely V = di ~ 
d/dx^. Formulas for the electromagnetic form factors Fi,2 
can be found in Appendix A. 

Having the Foldy-Wouthuysen transformation defined by 
the operator S (see Ref. and n = -^2(0)) 



S = — < /3 Q? • TT — 
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the new Hamiltonian is obtained via 
and takes the form 
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where by dots we denote the omitted higher-order terms, 
{X,Y] = XY + Y X,&nAa'j = [a\ crJ]/(2i). The Hamil- 
tonian H' is a generalization of the Foldy-Wouthuysen Hamil- 
tonian Hpw 1 24] to an arbitrary number of dimensions. The 
electromagnetic field in H' is the sum of the external Coulomb 
field, the external (constant) magnetic field, and a slowly vary- 
ing field of the radiation. For practical calculations, it is more 
convenient to have a Hamiltonian expressed in terms of the 
gauge-independent field strengths. To achieve this, we sep- 
arate out the Coulomb field and perform the Power-Zienau 
transformation of the Hamiltonian H' with the operator S' of 
the form 



S' = —e / dur ■ A{ur,t) . 
Jo 
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After neglecting irrelevant spin-independent terms, the trans- 
formed Hamiltonian becomes 
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Here, ^ p> - and B,k = V^B. H" is the gen- 
eralization of the Power-Zienau Hamiltonian Hpz Ii24il to an 
arbitrary number of dimensions. 

The Hamihonian H" includes most of the radiative cor- 
rections that are needed for our calculation, but not all of 
them. First, the higher-order terms with the anomalous mag- 
netic moment are omitted in H" . This contribution is more 
conveniently calculated with the exact Dirac-Coulomb wave 
functions, starting directly from the Hamiltonian (I10> . Fur- 
thermore, there is an additional correction that cannot be ac- 
counted for by the Fi and F2 form factors. It is represented 
by an effective local operator that is quadratic in the field 
strengths. This operator is derived separately by evaluating 
a low-energy limit of the electron scattering amplitude off the 
Coulomb and the magnetic fields. Details of this calculation 
are presented in Appendix B. The result is 



SH = — \2 g'^ B'^ ^^E^ 71 + a'^B'^ V^E^ f \ , (17) 
2 m'- ■' 

where B = const, E is an arbitrary electric field, and the func- 
tions T] and ^ are given by Eqs. (IB16> and (IB17t . respectively. 



III. ONE-LOOP SELF-ENERGY CORRECTION 

The dimensionally regularized NRQED approach formu- 
lated in the previous section will be first employed for a 
derivation of the self-energy correction to order {a/ir) {Z a)^ 
for the bound-electron g factor This derivation will serve us 
as a test of the new approach (as this result has been already 
obtained in our previous work fcoti ) and also as a basis for the 
two-loop calculation. 

As in 1I20I1 . we separate the one-loop self-energy correction 
up to the order of (a/n) {Z a)* into three parts. 



5« = 9'^' 
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where the first part is the the contribution due to the free- 
electron form factors Fi and F2, the second part is the con- 
tribution induced by the additional Hamiltonian and the 
third part is the contribution coming from low-energy photons, 
i.e. a Bethe-logarithm type contribution. 

We start with the form-factor part g'-^\ The anoma- 
lous magnetic moment ^2(0) in the modified Dirac-Coulomb 



Hamiltonian (I10> leads to the following energy shift linear in 
the magnetic field {d = 3) 
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where A = {B x f )/2, and we denote the one-loop com- 
ponents of the form factors by the corresponding superscript. 
The corresponding correction to the g factor is 
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In obtaining this result, we used the closed-form expression 
1I25I I26I1 for the component of the Dirac wave function per- 
turbed by the magnetic interaction, which has the same rela- 
tivistic angular momentum, as the reference state. 

For the remaining part of the form-factor contribution, we 
employ the transformed Hamiltonian M6\ . The last term of 
this Hamiltonian 



[F'l (0) + 2 F'2(0)] a'^ V^' [4 tt Z a 5\r)] B'^ r 
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gives rise to a contribution 

9?i - - (0) + 2 Ff^ (0)] (4 TT Z a <5'^(r)) . (22) 
The second-order correction to the energy 
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The total form-factor contribution is 
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(27) 



The second part of Eq. (I18> . denoted by gj^^ is a high- 
energy correction that is not accounted for by the form factors. 
It is given by the effective Hamihonian ( I17> . with E being the 
electric Coulomb field. The corresponding correction to the g 
factor is 
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where rj^^^ and are the one-loop components of the coef- 
ficient functions given below in Eqs. ( IB16> and ( IB17t . 

The third part of Eq. (I18> is a low-energy contribution that 
can be considered as a correction to the Bethe logarithm due 
to the interaction with an external magnetic field. Let us first 
derive the Bethe-logarithm correction to the hydrogen Lamb 
shift within the dimensional regularization. The correction to 
the energy is 
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Here, 6ljl = 5^^ — k^ k^ /k'^ is the transverse delta function, and 
k = \k\. After performing the integration over /c and dropping 
a common overall factor of (4 tt)^ r(l + e), SEl becomes 

a 1 (AnZaS'^ir)) 



SEl — 
TT be 



H_lln(Za)2]-llnfco 



where the Bethe logarithm In ko is given by 



hi ko = 



p{Ho - Eo) In 



2{Ho-Eo) 



{p{Ho-Eo)p) 



(30) 



(31) 



We now consider all corrections to 5El due to the presence 
of the external magnetic field. The first one is induced by the 
correction to the Hamiltonian [the fifth term on the right-hand 
side of Eq. ( I16H 
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The corresponding energy shift is given by 
SaE - ■ '^'^ 



Ho 
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where by Sa{- ■ ■) we denote the first-order correction to the 
matrix element induced by the perturbing Hamiltonian SaH. 
This matrix element is calculated using the scaling properties 
of the Schrodinger Hamiltonian given by Eq. (|9}, and the cor- 
responding correction to the g factor is found to be 
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The second correction to the interaction with the magnetic 
field is [sixth term in Eq. ( I16H 
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where the last part of the equation holds only for S-states. The 
corresponding contribution to the g factor is 
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The third correction is due to the coupling with the radiation 
field, [seventh term in Eq. ( I16H . 
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Here, the last expression is obtained by c?-dimensional angular 
averaging. The corresponding energy shift is written as 
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The contribution to the g factor is 
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The fourth contribution involves both the correction to the 
coupling with the radiation field and the interaction with the 
magnetic field, [the fourth and the ninth term of Eq. (I16H. 



and is of the form 
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The corresponding correction to the g factor is 

TT O c 



pling with the radiation field and the same interaction with the 
magnetic field, 
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The fifth contribution is due to another correction to the cou- 
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and is of the form 
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The corresponding correction to the g factor is 
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The sixth and the last contribution is due to the spin-orbit interaction and the interaction to the magnetic field, 
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This correction involves a more complicated matrix element with three propagators. 
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The corresponding correction to the g factor is 
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where In fcs is implicitly defined by the relation 
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which holds in the limit of large e. 

Finally, the total Bethe-logarithm type contribution to the g 
factor is a sum of calculated terms 

9f'=9'^l+9'^^+9'i^+9'^^.+9f^+9^l (50) 

The complete one-loop self-energy correction to the bound- 
electron g factor is then 
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in full agreement with the former result in Eq. (12) of 
Ref. 



IV. TWO-LOOP CONTRIBUTION 

The derivation of the two-loop corrections to the bound- 
electron g factor is performed in full analogy to the one- 
loop calculations. The total two-loop correction of the order 
{a/ IT? (Z oi? can be separated into four parts. 
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The first part is a form-factor contribution. The sec- 
ond part g\ is an additional high-energy contribution not ac- 
counted for by the form factors. The third part arises from a 
contribution in which one of the two virtual photons is of low 
energy. The second photon effectively modifies the vertex, 
which can be accounted for by the anomalous magnetic mo- 

(2) 

ment. The fourth contribution g\ involves the closed fermion 
loops and is called the vacuum polarization part. 

We start with the form-factor contribution. The two-loop 
anomalous magnetic moment correction is obtained from the 
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corresponding one-loop contribution, Eq. ( 1201 
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The one-loop anomalous magnetic moment in the modified 
Dirac-Coulomb Hamiltonian ( I10> leads to the following en- 
ergy shift. 
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The corresponding correction to the g factor is 
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The other contributions due to the two-loop form factors are 
immediately obtained from the corresponding one-loop ex- 
pressions in Eqs. ( I22> . ( I24> and ( I26> 

gfj = - \F'f (0) + 2 (0)] (4 TT Z a ^J'^Cr)) , (56) 

gf^ =2 (3 - 8 £) ^^'1'^ (0) (4 ^ Z a J'^(r)) , (57) 

gfi = - (1 - 4 £) F'l'^ (0) (4 ^ Z « <5'^(r)) . (58) 

The second-order corrections involving the slope of the one- 
loop form factors and the one-loop anomalous magnetic mo- 
ment vanish. It becomes clear if we notice that the coupling 
of the anomalous magnetic moment to the magnetic field, as 
obtained from the Hamiltonian jl6> . is 
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induced by the effective Hamiltonian 5H in Eq. ( I17> . with E 
being the electric Coulomb field. The corresponding correc- 
tion to the g factor is 
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The third part of Eq. i52\ g^^ is obtained from the formu- 
las for the one-loop Bethe-logarithm corrections. The overall 
coefficients in these formulas are modified by the presence of 
the anomalous magnetic moment n, in accordance with the 
corresponding terms in the effective Hamiltonian (16). The 
resulting corrections to the Hamiltonian describing the inter- 
action with the magnetic field are given by (for S-states) 
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The resulting two-loop corrections to the g factor are 
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The total two-loop Bethe-logarithm contribution is 
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All other possible two-loop corrections, which involve one- 
loop form-factors, are of higher order in the Z a expansion. 
Therefore, the total form-factor contribution is given by the 
sum 
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The second part of Eq. (|52} is a high-energy correction that 
is not accounted for by the form factors. This contribution is 
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The last part of Eq. i52\ g^ involves the vacuum- 
polarization correction. The contribution of the diagrams 
with the closed fermion loop on the self-energy photon is ac- 
counted for by the corresponding parts of the electromagnetic 
form factors Fi, F2 and 7/, ^. The two-loop vacuum polar- 
ization correction can be obtained from the correction due to 
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The corresponding contribution to the g factor is 
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The mixed self-energy and vacuum-polarization correction 
can be obtained in a similar way by the replacement 
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The corresponding contribution to the g factor is 
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The total vacuum-polarization contribution beyond the one ac- 
counted for by the form factors and 77, ^ is 
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Finally, the complete two-loop correction to the bound- 
electron g factor is given by the sum of four parts in Eq. \52\ . 
which yields 



7(2) 



a 

TT 

379 
"90" 



iZa) 



TT^ In 2 



28 
"9" 
379 
"6O" 



\n[{Za)-^ 
C(3) + - 



258917 
' 19440 
985 



In fco 



In 



113 
810 
5 



1728 144 



A ln2- — C(3) 
24 16 ' 



(71) 



The numerical values for In /cg and In for the first 7 5 states 
are 

In /coils') = 2.984 128 556, InfegllS) = 3.272 806 545, 

(72a) 

lnfco(25) = 2.811 769 893, lnfc3(2S') = 3.546 018 666, 

(72b) 

lnfco(35) = 2.767 663 612, lnfc3(3S') = 3.881960 979, 

(72c) 

lnfco(45) = 2.749 811 840, lnfc3(4S') = 4.178 190 961, 

(72d) 

lnfco(5S') = 2.740 823 727, lnfc3(5S') = 4.433 243 558, 

(72e) 

lnfco(65) = 2.735 664 206, lnfc3(6S') = 4.654 608 237, 

(72f) 

lnfco(7S) = 2.732 429 129, lnfc3(7S') = 4.849 173 615, 

(72g) 

The total numerical value of the nonlogarithmic term in 
Eq. (|7D for the 15 state is -16.436 842. All terms involving 
the closed fermion loop contribute —3.278 177 to this result, 
with the dominant contribution originating from the two-loop 

... . (2) 

vacuum-polarization correction g\J^ . 

V. RESULTS AND DISCUSSION 

In Table U we collect all contributions available for the IS* 
bound-electron g factor in three specific hydrogenlike ions 



which are important from an experimental point of view. For 
two of them, carbon and ojcygen, accurate experimental re- 
sults are presently available P.'s'], whereas the experiment on 
calcium is planned for the future |27]. 

The errors of the point-nucleus Dirac value and of the free 
part of the one-loop QED correction indicated in the table 
originate from the uncertainty of the fine-structure constant, 
a-i = 137.035 99911(46)121. The finite-nuclear-size cor- 
rection was re-evaluated in this work using the most recent 
values for the root-mean-square (rms) nuclear radii |28]. The 
error ascribed to this correction originates both from the un- 
certainty of the rms radius and from the estimated model de- 
pendence for the nuclear-charge distribution. 

The one-loop QED correction up to the order of {Z a)^ is 
given by the sum of the self-energy part [Eq. j51H and the 
vacuum-polarization part i29ll . 

,>'>=^(Z„,.(-H). (73) 

The higher-order one-loop self-energy correction was inferred 
from the results of the all-order numerical calculation 1 10, llil. 
For carbon and oxygen, the results presented in the table were 
obtained in Ref. |20] by an extrapolation of the numerical re- 
sults 1 10] for Z > 8, after subtracting the known terms of 
the Za expansion. The one-loop vacuum-polarization correc- 
tion consists of two parts, the electric-loop contribution that 
is due to the vacuum-polarization insertion into the electron 
line and the magnetic-loop contribution, which corresponds 
to the insertion of the vacuum-polarization loop into the in- 
teraction with the external magnetic field. The values for the 
higher-order electric-loop contribution presented in the table 
were inferred from the all-order numerical results of Ref. L30tl . 
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TABLE I: Individual contributions to the Is bound-electron g factor. The abbreviations used are as follows: "h.o." stands 
for a higher-order contribution, "SE" - for the self-energy correction, "VP-EL" - for the electric-loop vacuum-polarization 
correction, "VP-ML" - for the magnetic-loop vacuum-polarization correction, "TW" indicates the results obtained in this work. 
^^2^1/2 jjjg root-mean-square nuclear charge radius. 







12q5 + 


16q7+ 


"Ca 


Ref. 






2.4703 (22) 


2.7013(55) 


3.4764(10) 


[2^ 


Dirac value (point 


nucleus) 


1.998 721354 39(1) 


1.997 726 003 06 (2) 


1.985 723 203 7(1) 


Finite nuclear size 




0.000 000 000 41 


0.000 000 00155 (1) 


0.000 000 113 (1) 




1-loop QED 


(Za)° 


0.002 32281947(1) 


0.002 32281947(1) 


0.002 322 819 5 






(Zaf 


0.000 000 742 16 


0.000 001319 40 


0.000 008 246 2 


[3ri 






0.000 000 093 42 


0.000 000 240 07 


0.000 002 510 6 


[2ffl 




h.o., SE 


0.000 000 008 28 


0.000 000 034 43 (1) 


0.000 003 107 7(2) 


[iffi+[2ai 




ll.o., VP-EL 


0.000 000 000 56 


0.000 000 002 24 


0.000 000 172 7 


[30] 




li.o., VP-ML 


0.000 000 000 04 


0.000 000 00016 


0.000 000 014 6 


[13J 


>2-loop QED 


(Zaf 


-0.000 003 515 10 


-0.000 003 515 10 


-0.000 003 515 1 






(Zaf 


-0.000 000 001 12 


-0.000 000 002 00 


-0.000 000 012 5 


[3ri 




(Za)^ 


0.000 000 000 06 


0.000 000 000 08 


-0.000 000 010 9 


TW 




h.o. 


0.000 000 000 00(3) 


0.000 000 000 00 (11) 


0.000 000 000 (100) 




Recoil 


m/M 


0.000 000 087 70 


0.000 000 117 07 


0.000 000 297 3 


[161 




h.o. 


-0.000 000 000 08 


-0.000 000 00010 


-0.000 000 000 3 


[151 


Total 




2.001041 590 18(3) 


2.000 047 020 32 (11) 


1.988 056 946 6 (100) 



whereas the magnetic-loop contribution was taken from the 
recent evaluation ifisll . 

The {Z af and {Z a)^ parts of the two- and more-loop 
QED correction comprise the two-, three-, and four-loop con- 
tributions to the free-electron q factor, multiplied by a kine- 
matic factor of the electron IBlll . The {Z a)^ part of the two- 
loop QED contribution was derived in the present work. The 
uncertainty due to higher-order two-loop contributions was es- 
timated as 



(2) 

5h.o. 



(2) 



[{Zaf 



gW[{Zaf 



(74) 



(n) 

where ^ is the n-loop higher-order QED contribution and 
is the n-loop {Z af QED contribution. 

The nuclear recoil correction to first order in the mass ratio 
m/M but to all order in Z a was calculated in Refs. iHQ^l- 
The leading recoil corrections to order [ra/MY and am/M 
were derived in Refs. |32, 33] for a nuclear spin 7=1/2 and 
recently generalized for an arbitrary nuclear spin in Ref. fl^ . 

Based on the data presented in Table|I] we conclude that our 
evaluation of the one- and two-loop QED corrections to order 
{Z aY improves the accuracy of the theoretical prediction for 
carbon by an order of magnitude, as compared to the previous 
compilation fToll . The resulting QED contribution to order 
[Z aY turns out to be rather small for carbon and oxygen, 
as a result of a cancellation between the logarithmic and the 
nonlogarithmic parts of this correction [see Eq. j71H. For cal- 
cium, to the contrary, the numerical contribution of the two- 
loop [Z aY correction is large and of the same order as the 
[Z oY correction. This indicates that the perturbative Z a- 
expansion approach is no longer effective in this region of Z, 
and a direct all-order numerical evaluation would be highly 
desirable. 

It is remarkable that among different sources of the theoret- 
ical uncertainty for calcium, the error due to the higher-order 
two-loop QED correction is by far the dominant one. This 
means that, if the prospective experimental investigation of 



the bound-electron g factor in calcium is performed on the 
same level of accuracy as for carbon, namely 10~^, a com- 
parison of the theoretical and experimental results would al- 
low one to identify the contribution of the non-perturbative (in 
Z a) two-loop QED effects with a 10% accuracy. 

The comparison of the theoretical and experimental results 
for the IS* bound-electron g factor in carbon and oxygen 
yields the presently most accurate method for determination 
of the electron mass 1 2]. Based on the theoretical g factor val- 
ues presented in Table HJ we obtain the following values for 
the electron mass derived from the experiments on carbon |^ 
and oxygen (in atomic mass units): 



m(i607+) 



0.000 548 579 909 32(29), (75) 
0.000 548 579 909 60(41). (76) 



The uncertainty of these results originates from the experi- 
mental value for the ratio of the electronic Larmor precession 
frequency and the cyclotron frequency of the ion in the trap; 
the uncertainty due to the theoretical prediction is more than 
by an order of magnitude smaller and thus negligible. 
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APPENDIX A: ELECTROMAGNETIC FORM FACTORS 

We consider the form factors defined by 

-> = + ^F2(g2) [l^^] , (Al) 

where q is the outgoing photon momentum. The form factors 
are expanded in a up to second order, 

= l + Fi'\q^) + F['\q^), 
F,{q') = F^'\q') + F^'\q% (A2) 



where the superscript corresponds to the loop order, i.e. to the 
power of a. The results for the form factors expanded into 
powers of up to read (in Z) = 4 — 2e): 



F['\q-) 



F^iq^) 



1 



:t + 2e- 



1' 



1099 77 



11 


1 


5 




240 ~ 


40e ~ 


48^ 


•(- 


11 ^ 






V60 


^ 120^ 


)1 


7 



1296 
491 5 



144 



C(2) 



47 
576 



1440 
119 
"36" ^ 



C(2 



24 
2C(2) 



VP 



VP 

1721 
12960 ^ 72£2 



1 75 

3C(2) ln2-— C(2) 



1 



1 

487 



VP 



|-3C(2) In2+^C(2) 



- - -C(2) 
216 8^^ ' 



VP 



11 

^10 
3 

+ 4 
61 
40 



C(2) In 2 
C(3) 
C(2) 



23 



14731 
28800 



C(2) - §C(3) 



40 



C(3) 



533 
1080 



-C(2) 
10^^ ^ 



1637 



19 



VP 



5040 720 1 



— C(2) In 2 
14^^ ^ 



689 
1050 



C(2) 



(A3a) 
(A3b) 
(A3c) 



(A3d) 



The subscript VP denotes the contribution to the two-loop 
form factors which involves a closed fermion loop. 



APPENDIX B: THE LOW-ENERGY LIMIT OF THE 
SCATTERING AMPLITUDE 

In this section we describe the evaluation of the low-energy 
limit of the spin-dependent part of the scattering amplitude 
that gives rise to the effective Hamiltonian illl . The scatter- 
ing amplitude under consideration is schematically depicted 
in Fig.n where the leftmost graph is the "tree" diagram and 
the remaining graphs represent the tree diagram "dressed" by 
a self-energy photon. The two-loop diagrams are not shown 
explicitly; they can be obtained from the one-photon ones in 
a standard way. Each graph contains two interactions with 
the external field, one of which is the interaction with the ho- 
mogeneous magnetic field (a 7' vertex) and the other is the 
interaction with the Coulomb field of the nucleus (a 7*^ ver- 
tex). From the one- and two-loop scattering amplitudes we 
additionally subtract a tree amplitude with the vertices modi- 
fied by the electromagnetic form factors Fi and F2. This pro- 



cedure removes the part that is already accounted for by the 
Hamiltonian il6\ and leads to a simple polynomial expression 
for the resulting amplitude. 

FIG. 1 : Feynman diagrams representing the scattering amplitude of 
a free electron on both the Coulomb and the magnetic field, at the 
tree and the one-loop level. 

In order to extract the spin dependent part of the scattering 
amplitude, we construct the projection operator. Let us first 
consider a general non-relativistic operator Q, 

Q = Q° + Q'a\ (Bl) 

The spin-dependent part of Q can be retrieved by the follow- 
ing projection operator: 

Q'^^TiiQa']. (B2) 
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In d dimensions, the nonrelativistic expansion of the Hamil- 
tonian involves ct*-' — [cr* ,(T-']/(2i). The extension of the 
spin-projection operator to an arbitrary number of dimensions 
is 



(B3) 



with Q ~ Q'^ cr'^ . We assume here the following properties 
of the trace to hold: 



Tr[CT'^] ^ 0, 
Tr[I] = 2, 



(B4) 



We now consider the operator Q sandwiched between the 
positive-energy solutions of the free Dirac equation normal- 
ized by = 1. The following identity holds, 

u{pi, Sf) Q u{pu Si) = Tr[(3 s{) u{pf, Sf)] . (B5) 

Since our aim is to calculate the low-energy limit of the am- 
plitude only, we can use an approximate form for u{p, s), 

( \ 
u{p, s) « , (B6) 

V -p^iis) J 

where is a nonrelativistic spinor Using a replacement that 
extracts the spin dependence 

0(si)0+(sf)-^^, (B7) 
the projection operator becomes (in units m = 1) 

u{pi,Si)u{pf,Si) i 



^a-piff'^, 



- \ cr*-' t? • Pi 



Therefore, 



+ 1) S*^' {^f + 1) , 



(B8) 



16 



Tr[(|<f + l)Q(j^i + l)S]^^] . (B9) 



We now turn to the scattering amplitude of the free electron 
on the Coulomb and magnetic fields. The spin-dependent part 
of this amplitude is written as 

Q^Q^'''PeAo{qi)eA^{q2)a.p, (BIO) 

where qi and q2 denote the exchange momenta. The ampli- 
tude corresponding to the tree diagram in Fig.^is given by 



1 



■Tr 



16 



+ 1) 7" ■ 
1 



1 



1 



r 



7°(^^i + l)S'''' 



(BID 



M ii - 1 

where the momenta pi, pi are on the mass shell, and the ex- 
change momenta are spatial, <Zi = <Z2 — 0- 

As an example of one-photon contributions, we give an ex- 
pression for the rightmost diagram in Fig.Q] 



d^k 1 1 



{2tt)D fc2 16 
1 



Tr 



1 



■r 



The other one- and two-loop contributions are obtained in the 
analogous way. From the resulting amplitude we subtract the 
tree amplitude with vertices 7° replaced by F", 



1 



■Tr 



+ i)rO(qi) 



16 

H^i + l)T^{q2)- 



■T^{q2){A + ^)^ 



- 1 



T\q,){^, + l)ll 



up 



(B13) 



where F" is defined in Eq. (lAH . The final expression for the 
total amplitude Q^^^f is obtained by the expansion in small 
momenta pi , pf and the subsequent integration over the loop 
momenta. The result for Q^-'^p can be written in the form 



(B14) 



where the functions T^'^p and Q^^'^P are orthogonal to ql^ (due 
to the gauge invariance) and antisymmetric in v, p. Their ex- 
plicit expressions are 



TP^P = qP {qPq!^ - q^qP) + q, ■ q^ {gPPql ~ gP^q^) , 
gP'-P = qligPPq^-gP'q'^) ■ (B15) 

The results for the coefficient functions rj and ^ read 



a 2 



^ 47r3e^U7r/ 
283 169 



/2528 169 , 

■TT" 



V 81 



54 
2 



VP 



— + TT^ - A TT^ In 2 + - C(3) - - 

10 120 15 5 3e 



(B16) 



^ 3£ j ^ V4^ 



'2674 91 



V 81 27 
102 



VP 

4 



152 319 2 68 2, „ 

\ TT^ TT^ In 2 + C 3) H 

15 45 5 5 ^ 3e 



(B17) 



where the subscript VP denotes the contribution involving 
a closed fermion loop. The effective local operator Q in 
Eq. ( IB10> becomes 



Q= -ViTP'^P + ^gP-'P]eA^{q^)eAM^'^P 

^ \2a'^ B''''V^E''r] + a'^B''^V''E''(\ , (B18) 

2 m'- ■' 



which corresponds to the effective Hamiltonian in Eq. ( I17> . 



X 74:^i + 1) s-^" 



+ symmetrization . (B12) 
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